Compact Stars in the QCD Phase Diagram III (CSQCD III) 
December 12-15, 2012, Guarujd, SP, Brazil 
http : //www . astro . iag . usp . br/~f oton/CSQCD3 



Explosive Combustion of a Neutron Star into a 
Quark Star: the non premixed scenario 



. . . Rachid Ouyed ^, Brian Niebergal ^ and Prashanth Jaikumar ^ 

^^ ■ ^Department of Physics and Astronomy, University of Calgary, Calgary, Alberta, 

(^ • T2N 1N4-, Canada; Email: rouyedSucalgary . ca 

^ ! ^Department of Physics and Astronomy, California State University, Long Beach 1250 

^ Bellflower Blvd., Long Beach, CA 90840-9505, USA 

<: 
o: 

! Abstract 

ftJ ' We review aspects of the hydrodynamical combustion of nuclear matter to 

hH . strange quark matter in a neutron star. Numerical studies on non-premixed 

: \ combustion that consistently include hydrodynamical flows in a reactive-diffusive 

o ' setup show that in ID, the conversion (burning) front moves at sub-sonic speeds 

JL . and stops short of converting the entire star to SQM, essentially due to advec- 

i-{ [ tive forces. However, in the process, we also find that neutrino cooling of the 

c/3 ' interface causes it to wrinkle, laying a platform for a deflagrative-to-detonative 

I transition (DDT). We outline progress on improvements in the burning code 

! (Burn-UD: | http://quarknova.ucalgary.ca/software/Burn-UD/[ ) that will ulti- 

^ I mately reveal the mechanism that can explode the outermost layers of even a 

OO ■ dense compact object like a neutron star. 

^' 

o: 

.• ' 1 Introduction 

^ . 

en ' The traditional strange-quark matter (SQM) hypothesis is the assertion that bulk 

matter (large baryon number) comprised of deconfined up, down, & strange quarks 

.^ . (or more conveniently written throughout the rest of this paper as u, d, & s) at zero 

^ ! pressure is more energetically favorable than the most stable atomic nuclei [H l2l [3] . 

j^ I Thus, once formed, SQM may exist forever. Nuclei are then only a metastable state. 

This hypothesis can be described using the MIT bag model - a simple model that 

nevertheless has been reasonably successful in describing thermodynamical properties 

of quark matter (see |1] and references therein) . 

In this paper, we review our numerical investigation of the issue of combustion 
of pure neutron matter to u,d,s matter. For SQM burning, far from turbulence- 
flame interactions, even the laminar flow has not been analyzed in sufficient detail. 
An improved prescription of the burning front in the laminar flow approximation 
was given in [5l|6]. Our simulations use hydrodynamics, taking into account binding 
energy release and neutrino emission across the burning front - going beyond previous 



treatments of the problem. Our studies already find speeds as high as ~ c/100, where 
c denotes speed of light. We find indications that a de-leptonizaton instability (a new 
type of instability unique to this type of conversion) and probably turbulent effects 
may well decide the fate of the conversion (deflagration or detonation). 

This problem is interesting for three main reasons: (i) in Type la supernovae, 
mult i- dimensional studies of small scale dynamics of flame burning including tur- 
bulence provide support for a pathway to the distributed regime, which can be a 
platform for detonation of the white dwarf [71 [HI E]. A delayed detonation model for 
Type la explosions has also been discussed in the context of material unburnt by the 
initial deflagration ^U\. Furthermore, deflagration-to-detonation (DDT) transition 
through turbulent flame burning in laboratory experiments on combustible gas mix- 
tures has been numerically studied [H] - all this naturally leads to an investigation of 
the existence of similar effects in SQM burning (ii) recent work on burning neutron 
matter to SQM p[2] uses simple scaling arguments to show that the above pathway to 
explosive detonation may occur for SQM burning in a neutron star and (iii) explosive 
conversion of a neutron star to a strange quark star is astrophysically relevant as 
discussed in this paper. 

Following this very brief introduction, the rest of the paper is structured as follows: 
In Section 2, we present results from numerical simulations of the combustion from 
hadronic to quark matter and provide a recipe on how to set up such simulations. 
This is essentially an investigation on how hadronic matter would transition to u,d,s- 
quark matter through the conversion of hadrons to M,(i-quarks, then u,(i-quarks to 
s-quarks. We also present an analysis of the similarities between this transition and 
a combustion process, specifically a combustion of the non-premixed type, a distinc- 
tion that is critically importanto. A discussion of the results and future outlook for 
numerical studies of the SQM burning are presented in Section 3 before we conclude 
in Section 4. 

2 Simulations of The Hadronic-to-Quark-Matter 
Phase Transition 

Although asymptotic freedom would appear to suggest otherwise, even at supra- 
nuclear densities (i.e, higher than that of a nucleus), nucleons cannot spontaneously 
dissolve into their constituent u & d-quarks as the overall energy per baryon would 
be higher. The presence of s-quarks adds an extra (flavor) degree of freedom, thus 
lowering the overall energy per baryon |T] ; resulting in the SQM phase of matter. Since 
the probability of converting many u- & (i-quarks into 5-quarks at the same time is 



^Previous literature on this topic has reported very different results on the transition speed and 
energy as a result of (incorrectly) assuming premixed combustion 



The Interface 



Q s-quark diffusion ^ Mean free path 
^^^^^^^^^^^^^^^' from scattering rate 




hadronic matter 



u,d,s-inatter 



O d-quark diffusion 



Figure 1: The interface where hadronic matter is one side and SQM on the other, with 
the interplay between reactions and diffusion governing the speed at the which the 
interface spreads. The presence of s-quarks activates the burning and so the process 
is governed by particle diffusion rather than heat diffusion. 



effectively zero, a conversion from hadronic to SQM can only proceed in practice if a 
substantial number of s-quarks are present. 

Consider the region inside a neutron star (NS) where the central density has 
reached that of nuclear deconfinement. Subsequently, a seed of SQM is instantiated 
at the star's center. It has been speculated that such a situation could occur during 
the core-collapse phase of a supernova, or, if the nucleation timescale is slower [13], in 
an older NS whose central density has increased due to spin-down [T3] . The details of 
the seeding mechanism are not considered here but there are several possibilities, for 
example clustering of lambdas, higher-order neutrino "sparking" reactions, or seeding 
from the outside [15] . 

In the boundary region surrounding the SQM seed, hadrons will then overlap with 
s-quarks, allowing the hadrons to dissolve into their constituent u- & d-quarks. Such 
a region will attempt to equilibrate chemically by producing more 5-quarks, which 
will diffuse into the hadronic regions, forming even more SQM. The situation is then 
an interface where hadronic matter is on one side and SQM on the other, with the 
interplay between reactions and diffusion governing the speed at which the interface 
spreads (see Figure [1]). This scenario is of a type similar to that of non-premixed 



combustion [16j, rather than premixed. The important difference is that the reactions 
in question are not activated above a certain temperature, as is often the case for 
combustion processes. Instead, the presence of 5-quarks activates the burning and 
so the process is governed by particle diffusion rather than heat diffusion. If one 
uses only the jump conditions, one is assuming that only fluid pressure differences 
matter - this is like setting fire to a barrel of oil, where everything is fiammable on 
contact. However, we take into account diffusion of strange quarks, starting from a 
non-premixed state, and this slows down the burning speed, but not by more than an 
order of magnitude surprisingly. Further, the neutrinos which change the temperature 
profile at the interface and alter the jump conditions quite a bit. Further issues 
regarding the nuances of the burning regime are discussed in § I2.3.2[ 

During the combustion an interface is created with cold (unburnt) fuel on one 
side and hot (burnt) ash on the other. Subsequent cooling of the ash will result 
in pressure gradients that induce fluid-like motion of the matter in and around the 
interface. To conserve energy, fluid velocities across the interface would then dif- 
fer, causing compression or rarefaction in the burning region which can enhance or 
quench the combustion. Thus, a numerical fiuid-dynamical treatment of the problem 
is appropriate, even essential, to determine the interface speed. 

This problem was addressed numerically in one dimension by [5]. We refer the 
interested reader to [6] where the problem is first formulated in a reactive-diffusive 
hydrodynamic setup (see § 4 in [6]). The importance of conserving energy and gener- 
ating entropy while re-arranging particles into different states is emphasized from the 
start. Details of the numerical procedure and relevant techniques are given in § 5 in 
[6] and instructions are included for obtaining and running Burn-UD. In § 6 in [6\, ba- 
sic results of Burn-UD such as the combustion front speed are verified with analytic 
solutions (jump conditions). The broader implications of the numerical and analytic 
results are described in § 7 in [6J , with applications to compact stars, relativistic heavy- 
ion collision experiments, and the early Universe. Finally, animations of the velocity 
and temperature profiles, as well as the reaction rate, neutrino emissivity and mean 
free paths can be found online at: "http://quarknova.ucalgary.ca/software/Burn-UD/|' 

The most notable results of the modern numerical approach are the rapid burning 
speeds (0.002-0.04 times the speed of light) and the quenching of the burning due to 
neutrino cooling, resulting in a new type of instability as discussed in §2.21 

2.1 Evolution of the chemical potentials 

The spatial profile of the different particle species is shown in Figure |2l Just after 
contamination by s-quarks, in zone B (cf. Fig. |2]), at small values oi fig, the reactions 
producing s-quarks are dominated by the (/i^ — /^s)^ factor in equation 7 in [5] (or eq. 
4.16 in [0]). Further behind the interface, s-quark production becomes increasingly 
dependent on the temperature term rather than the chemical potential difference. 
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Figure 2: A snapshot of the chemical potentials for different species [u,d,s quarks, 
and electrons) throughout the burning interface. Meta-stable u,d matter (fuel) is 
upstream of the interface (zone A). Upon contamination by diffusing s quarks (oxi- 
dizer), the non-leptonic reaction (eq. 4.15 in |6]; eq. 6 in [5j) burns the d quarks into 
s quarks (zone B). In zone C the leptonic reactions (eqns. 4.11-4.14 in |E]; eqns 2-5 in 
[5]) become the dominant means of producing s quarks, leading finally to equilibrated 
SQM. The relevant initial parameters are density /xinit = 300 MeV and electron frac- 
tion Xe = 0.04nTOT- Leaving hydrodynamics out of the simulations has little effect on 
these abundance profiles but enlarges the width of the interface by a factor of ~ 10. 



which tends to increase the burning speed. This importance of temperature on the 
burning process has also been realized in analytical studies (ITJ ITS] . 

The s-quark mass value does have a mild effect on the burning speed in the 
following way. The front speed, as we will show in §2.2[ is tempered by neutrino 
cooling, which in our simple treatment of neutrino transport, depends mainly on the 
neutrino mean free path due to scattering on electrons. The electron density drops as 
the fuel (M,(i-quark matter) is burnt. A larger value for the s-quark mass leaves behind 
a greater electron fraction, resulting in less neutrinos leaving the system. This implies 



that the interface will move faster, however, this effect is negligible compared to the 
theoretical uncertainties in the neutrino cooling rate. Thus, changing the s-quark 
mass, or the initial electron fraction, modifies the width and speed of the interface 
only by a few percent. 

The initial temperature also plays a minor role, provided the assumption of cold 
degenerate matter (T -C /x ~ 300 MeV) is maintained. Hence, an independence on 
all initial parameters, except density, is exhibited by the system. 

The bag constant was kept fixed even though the behavior of the interface depends 
upon its value. The justification is that the bag constant may have ill-described 
dependencies on other variables, for example the chemical potential or temperature 
[19]. A more self-consistent study then would be one that, rather than varying the 
bag constant, explores different equation-of-states. 

2.2 A New Type of Instability: wrinkling due to de-leptonization 

The inclusion of neutrino cooling in the simulations, through the production and 
transport of neutrinos (as described in §4.1.3 in |6]; see also |5] ), results in a reduced 
pressure for the burnt fluid. This gradient forces a decrease in fluid velocities both 
behind (f 2) and in front (f 1) of the interface, causing advection to oppose the progres- 
sion of the burning interface. Hence, as the neutrino cooling rate increases the burning 
speed decreases. An analytic treatment using the jump conditions (energy-momentum 
conservation) confirms this effect and also predicts that, for a given density there is 
a critical neutrino-cooling rate that quenches the burning [5]. 

Above the critical cooling rate, instead of moving backwards the interface comes to 
a halt, since as the interface stops diffusing into the fuel, reactions no longer proceed 
and neutrino production ceases. Thus, energy is no longer being removed from the 
burnt region, which causes the pressure gradients to weaken. The system then reaches 
a situation where diffusion and advection are in balance, and the interface is stopped. 

In other words, the velocity of the (backwards-moving) fluid ahead of the interface 
is proportional to the neutrino cooling, unless it becomes larger than the reactive- 
diffusive burning speedcl f rd < I'^il; at this point it loses its dependence on neutrino 
cooling and becomes constant. Only at lower densities (< 2 times nuclear saturation 
density) is the interface moving slow enough to allow the burnt material to cool suffi- 
ciently and halt. Otherwise the neutrino cooling rate would need to be unphysically 
large to halt the combustion. 

Consider now the interface in more than one dimension. As a few regions along the 
burning interface are halted due to cooling, unburnt material starts to flow backwards 
onto the interface, whereas regions not halted by cooling will be pushing unburnt ma- 
terial away from the interface. The burning interface is mostly opaque to neutrino 



"The reactive-diffusive burning speed, vjio, is the burning speed in the absence of any fluid 
velocities. 
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Figure 3: A snapshot during the simulation of the temperature (T ), reaction rate 
(-R), and neutrino emissivity, e^, throughout the burning interface. The temperature 
is shown with (sohd hne) and without (dashed hne) neutrino coohng effects, where 
the difference between the two is the variable C = T — Tcooied that serves as the 
measure of cooling. 



cooling, but not entirely. There is a small drop in temperature behind the inter- 
face (see Figured]). Since Pressure goes as T^ this is significant. The burnt areas 
end up underpressured, which halts the progression of the interface in that region. 
The neighbouring regions are still advancing, which creates severe wrinkling, which 
increases the burning speed. The result is a wrinkled interface, with a shearing be- 
tween the unburnt fluids (fuel) in halted and non-halted regions. Since a wrinkled 
interface increases the diffusion rate, the effect is an overall increase in the burning 
speed. Thus, the non-halted regions will tend to burn even faster, thereby increasing 
the wrinkling - leading to an instability. 

The wrinkling will be subject to stabilization - in the dimension along the interface 
diffusion causes concave regions to accelerate while convex regions are decelerated. 
Whether this stabilization is able to overcome the instability, or, if the instability 



grows unhindered is an open question requiring high-resolution multi-dimensional 
simulations. In any case there is a new type of instability which, due to deleptoniza- 
tion, can halt the combustion of u,d-qua.Tk matter. 

2.3 Combustion of a Neutron Star 

In the case of a NS burning into u,d,s-quaik matter, there is a radial density gradient 
in the star that the burning interface progresses through. As the interface burns 
outwards it reaches lower densities, making the interface move slower, until it reaches 
the critical point where the interface halts due to the above-mentioned deleptonization 
process. Note that, the critical point depends on a number of factors including the 
equation-of-state, neutrino cooling efficiency, and most importantly the dynamics of 
the combustion in multiple dimensions (eg. the wrinkling of the combustion front, 
sensitivity to instabilities, etc.). 

2.3.1 Low-Pressure Cavitation 

Eventually the entire interface will reach low density regions and halt. The core will 
be underpressured, and the outer material will be falling into it (being burnt as it 
does). The situation is then an under-pressured u,d, s-quark core, with the outer layers 
of the NS lying on top. The under-pressured core will eventually collapse - similar to 
the core collapse phase during a supernova - causing the outer lying material to fall 
onto the core, leading to a second explosioiu. This sequence of events is illustrated 
in Figure HI 

2.3.2 Combustion Regime, Turbulence, &: Jump Conditions 

As noted above, the burning regime is that of non-premixed combustion [16], rather 
than premixed. Although this distinction may seem trivial it has important conse- 
quences, such as invalidating assumptions made in the literature for turbulent burn- 
ing, specifically the fiamelet-regime approximation. We expect turbulence to have a 
significant impact on the wrinkling of the combustion interface. This will increase 
the burning speed, not only due to an increase in the surface area of the combustible 
materials, but also due to a significant increase in the diffusion. An important point 
to realize is the effect of temperature/cooling on the combustion to quark matter. Not 
only can cooling change the properties of the interface such as width and speed, but 
can quench the burning process entirely. Thus, it is essential to include temperature 
increase when solving the jump conditions, as was done in [5]. There are numerous 
examples in the literature where the jump conditions are solved and the temperature 
of the SQM phase was set to zero! On the issue of turbulence, we remark that Herzog 



^Thc first explosion being the supernova that created the NS. 




Figure 4: An illustration of the burning of a NS. Region A denotes the initial burning, 
where the interface is laminar and moving outwards at 10^'^ to 10~^ times the speed 
of light. In region B instabilities - such as the Rayleigh- Taylor instability - are likely 
to develop, causing the interface to wrinkle and proceed even faster. At this point 
either the interface will detonate and eject the outer layers of the original NS, or, if 
the wrinkling is not strong enough, the interface will be subject to the deleptonization 
instability (see §2.2p and halt (region C). The outer layers of the original NS will then 
fall onto the newly created u,d, s-quark core, leading to a core-collapse and subsequent 
core-bounce. In either case a highly-explosive event is realized. 

and Ropke [20] recently presented results on the burning of neutron matter to SQM 
including turbulent effects such as Rayleigh- Taylor instabilities. This is distinct from 
the neutrino-driven instability we mention above. Qualitatively, they too find that 
the burning front stops because the reactions are no longer exothermic. However, 
this ignores the possibility of having a transition to a detonative regime driven by the 
neutrino-driven instability. 
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2.3.3 Transition to Superconducting Phases 

In a manner similar to Cooper pairing of electrons, weakly interacting quarks at 
high density are expected to pair in order to minimize their free energy, resulting 
in a color-superconducting state (e.g. the Color-Flavor Locked phase; [21] )■ This 
superconducting phase transition is expected to occur in quark stars shortly after they 
are formed, since cooling processes are very efficient in lowering the star's temperature 
below the rather large critical temperature for color superconductivity T^ ~ 10 MeV. 
The effects of color superconductivity can be taken into account upon the formation of 
a superconducting interface, where on the superconducting side an extra term (oc A 
the pairing energy gap) would have to be included in the equation-of-state. The 
Burn-UD could then be used to study the dynamics of this interface. Care would 
need to be taken in order to ensure electrodynamic constraints across the interface 
are satisfied (further details are given in Part III in [6J). 

Furthermore, in one of these phases, the Color-Flavor Locked (CFL) phase, the 
photon emissivity dwarfs the neutrino emissivity [221 l23]. In this situation matter 
at a certain radius in the stacl will drop below Tc causing matter in this region to 
release energy - (photons and neutrinos) freed as a result of the decay and interaction 
of low-energy (scalar and vector) modes. One can obtain a very high number density 
of thermal photons (photon fireball) that can drive the outer layers of the NS like a 
piston, with important consequences for heavy-element nucleosynthesis [2] 



3 Discussion and Outlook 

• Nucleation: A more detailed picture of the pathway to the quark-hadron phase 
transition is required if we are to establish astrophysical implications of the tran- 
sition. Basically, one needs to know how long it takes for a seed of quark matter 
to nucleate and grow inside the NS before it can be analyzed as a bulk fiuid 
(as was implicitly assumed in the combustion studies). Knowing the nucleation 
time accurately is critical. For example, in several long GRBs, an association 
with a prior supernova event is established, with time delays between few days 
to several months (e.g. [251 126] and references therein). This could be indicat- 
ing the delayed onset of the quark-hadron phase transition, which can release 
sufficient energy to power the gamma-ray burst when neutron matter is burned 
to quark matter. Recently, this scenario has also been applied to detailed light- 
curve modelling of extended emission, afterglows and late-time X-ray fiaring 
detected in several long GRBs by the Swift sateUite [23- Thermal nucleation 
timescales are also relevant to the proto-NS phase after the supernova explodes. 



^It is not clear whether the inner or outer layers of the star will drop below the critical temperature 
first 
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when the appearance of a quark phase can delay the collapse to black holes 
until after deleptonization is complete. Updated calculations of the nucleation 
rate and timescale taking into account these modified expressions for the shear 
viscosity of neutron matter and kaon condensed matter. The earliest studies of 
nucleation used the framework of quantum nucleation with the WKB approxi- 
mation for tunneling, obtaining nucleation times in cold /3-stable and hyperonic 
matter that had large uncertainties [Tnuc is exponentially sensitive to the over- 
pressure). We are presently working to reduce some of these uncertainties using 
updated transport coefficients in dense matter. 

The Deflagration-Detonation- Transition (DDT): The question remains 
open as to whether the combustion will switch from a deflagration to a detona- 
tion. Given the rapid flame speed that we find {v ~ 0.01c) and the size of the NS 
i? ~ 10 km, the conversion timescale of order 10"'^ seconds is of the same order 
as for the development of well-known flame instabilities - the Landau-Darrehius 
and Rayleigh- Taylor instabilities. Therefore, we have good reasons to expect 
that the laminar flow analysis in ID is insufficient and turbulent flame burning 
must arise in 2D. Given this, it would seem likely to see a substantial increase in 
burning speed, but results from high-resolution mult i- dimensional simulations 
are needed to make any conclusions. However, we note that in either case a 
second explosion would be realized; either the detonation explodes the outer- 
layers of the parent NS, or, the interface halts (due to the newly discovered 
deleptonization instability) and the outer-layers fall back onto the SQM core 
resulting in a core-bounce explosion. The energy of the explosion is then due 
to either to chemical energy (the detonation), or to gravitational energy from 
core-collapse. 

Temperature effects on SQM burning: 

Another important discovery was the effect of temperature on the combustion 
to quark matter. Not only can cooling change the properties of the interface, 
such as width and speed, but can quench the burning process entirely. Thus, it 
is essential to include temperature increase when solving the jump conditions, as 
was shown in [3] (see also § 6 in [B]). It is important to keep in mind that while 
the Bag model (equation-of-state) has been extremely successful at explaining 
observations, both in experiment and in Nature, it is only considered to be 
a rough approximation. QCD allows for more sophisticated models, some of 
which may change the dynamics of SQM burning. The Burn-UD could readily 
be used to test these models but is left as an avenue for future work. 
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4 Conclusion 

We have performed ID numerical simulations of the burning of neutron matter to 
strange quark matter with consistent treatments of reactions, diffusion, and hydro- 
dynamics. Our two new results, viz., the generation of a rapid conversion front and 
its possible instability due to cooling effects by neutrino emission are exciting, but 
only suggestive at this point. Clearly, more sophisticated numerical work is required 
to verify our findings in a multi-dimensional setting, but there is widespread interest 
in modelling QCD phase transitions inside compact stars/core-collapse supernovae 
and the early universe. Such groups could gain insight from our results and perform 
additional simulations. 

The Burn-UD code is publicly available, written in a modular form allowing it 
to be easily adapted to other problems of flame burning and propagation in non- 
premixed fluids. The explosion from the transition of hadronic to u,d, s-quark matter 
is not yet proved from first principles, but there are hints from the ID study that 
such an explosion is possible. If confirmed, this would vahdate the Quark-Nova modeO 
[28l 129] . consequences of which include: i) nucleosynthesis of rare elements [21], ex- 
planations for; ii) gamma-ray bursts [301 123] ; iii) ultra-high energy cosmic rays [31] ; 
iv) soft-gamma repeaters, anomalous x-ray pulsars, x-ray dim NSs [321 ESI EH ESI ES] , 
and ; v) super- luminous supernova [23 ESI ESj- The vahdation of the Quark- Nova 
would also have implications to Cosmology [101 HH |12l US]. 

R. O. is grateful to the organizers of the CSQCDIII for an inspiring meeting. This 
research is supported by a grant from the National Research and Engineering Council 
of Canada (NSERC). 
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